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We present theoretical analysis and numerical studies of the quantized vortices in a rotating 
Bose-Einstein condensate with spatiotemporally modulated interaction in harmonic and anharmonic 
potentials, respectively. The exact quantized vortex and giant vortex solutions are constructed ex- 
plicitly by similarity transformation. Their stability behavior has been examined by numerical 
simulation, which shows that a new series of stable vortex states (defined by radial and angular 
quantum numbers) can be supported by the spatiotemporally modulated interaction in this system. 
We find that there exist stable quantized vortices with large topological charges in repulsive con- 
densates with spatiotemporally modulated interaction. We also give an experimental protocol to 
observe these vortex states in future experiments. 
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INTRODUCTION 

The investigation of rotating gases or liquids is a cen- 
tral issue in the theory of superfluidity [lH4[ since rota- 
tion can lead to the formation of quantized vortices which 
order into a vortex array, in close analogy with the be- 
haviour of superfluid helium. Under conditions of rapid 
rotation, when the vortex density becomes large, atomic 
Bose gases offer the possibility to explore the physics of 
quantized vortices in novel parameter regimes. During 
recent years, there have been advances in experimental 
discoveries [5|4Z| of rotating ultra-cold atomic Bose gases, 
and these developments have been reviewed in [8j . 

Theoretical studies mainly make use of the mean-field 
Gross-Pitaevskii (GP) equation to describe the main fea- 
tures of the vortex states 0, Hoj |, and several predic- 
tions [To[ have been shown to agree with experiments 
1 ll | . Some of the important studies were concerned with 
the equilibrium properties of a single vortex, including 
its structures and dynamics [To|, the critical frequency 
and the nonlinear dynamics of vortex lattice formation 
[l2| . A multi-quantum vortex is typically dynamically 
unstable in harmonically trapped BEC predicted by sev- 
eral theoretical studies [13l4l5|. The splitting instability 
in case of multi-quantum vortices shows that the vortex 
will split into single quantum vortices even in the absence 
of dissipation due to the peculiar feature of nonlinear dy- 
namics Oil- However, in the presence of a plug potential 

Q, var- 
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[17j or an anharmonic trapping potential 
ious studies have addressed different means to stabilize 
multi-quantum vortices in rotating BEC. For example, 
when the confining potential is steeper than harmonic 
potential in the plane perpendicular to the axis of rota- 
tion, multi-quantum vortices are energetically favorable 
if the interaction is weak enough. For stronger inter- 
actions, the multiply quantized vortices break up into 
arrays of several vortices. In addition, interestingly, sta- 



ble multi-quantum vortices have also been found to exist 
in two-component BEC [20|, which can be adjusted near 
Feshbach resonance through spatial inhomogeneous ex- 
ternal magnetic field i.e. B = B(x). 

Mathematically, the GP equation, to be written ex plic - 



. explic- 



itly, is an equation of nonlinear Schrodinger type 
This equation has been studied extensively both in the 
physical and mathematical literature, since they provide 
a universal model for a study of the dynamics of the en- 
velope waves. One of the distinctive features of the equa- 
tion as it appears in BEC problems is the presence of an 
external trapping potential, which essentially affects the 
elementary excitation spectrum. Most properties of the 
BEC are significantly affected by the interatomic interac- 
tion, which can be characterized by the s-wave scattering 
length [28[- Recent experiments have demonstrated that 
both amplitude and sign of the scattering length can be 
modified by utilizing the Feshbach resonance [29|. This 
technique provides a very promising method for the ma- 
nipulation of atomic matter waves and the nonlinear ex- 
citations in BEC by tuning the interatomic interaction. 
By using this technique, one can study atomic matter 
waves and the nonlinear excitations in BEC for the case 
of the GP equations with the time- and space-dependent 
nonlinearity coefficients [22|-l25|. 

Motivated by stabilizing multiple vortex states and 
understanding the behavior of nonlinear excitation in 
physical systems, we perform theoretical analysis and 
numerical studies of the quantized vortices in a rotat- 
ing BEC with spatiotemporally modulated interaction in 
harmonic and anharmonic potentials, respectively. Com- 
pared with the former work on quantized vortices, we 
find that a new series of exact single and multiple vortex 
states (defined by radial and angular quantum numbers) 
can be supported by the spatiotemporally modulated in- 
teraction in a rotating BEC. In particular, our results 
have provided a very promising method for stabilizing 
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the vortex having very large topological charge S > 2, 
which has been conjectured unstable [30| by tuning the 
external potential and nonlinear interaction simultane- 
ously in time. 



THE THEORETICAL MODEL AND EXACT 
VORTEX SOLUTIONS 

At zero temperature, the quantum and thermal fluctu- 
ations are negligible so that a BEC trapped in an external 
potential can be described by a "macroscopic wave func- 
tion" ^f(r,i) obeying the GP equation. In the rotating 
frame with rotating frequency Qq around the z-axis, the 
GP equation in cylindrical coordinate reads [? ] 

dfy fr 2 h 2 d 2 o 

where V 2 = d 2 /dx 2 + d 2 /dy 2 = d 2 /dr 2 + 1/r x d/dr + 
l/r 2 d 2 /d0 2 with r 2 = x 2 + y 2 , m is the atom mass, 
is the azimuthal angel, the wave function is normal- 
ized by the total particle number N = Jdr|^| 2 , V ex t is 
an external trapping potential, and G = 47rh 2 a(r,t)/m 
represents the strength of interatomic interaction char- 
acterized by the s-wave scattering length a(r, t), which 
can be adjusted experimentally by an inhomogeneous ex- 
ternal magnetic field B = B(x,y,t) in the vicinity of 
a Feshbach resonance [29]. The trapping potential can 
be assumed to be V ext = m(uo 2 r 2 + u 2 z 2 )/2, where uo r 
and uo z are the confinement frequencies in the radial and 
axial directions, respectively, and in particular, the ra- 
dial confinement frequency uo r is assumed to be time- 
dependent as in (13, [3l|. In the following, we consider 
the atoms in the \F = l^mp = 1) hyperfine state of 7 Li 
and \F = l^mp = 1) hyperfine state of 87 Rb trapped in 
a very thin disc-shaped potential, i.e., the trapping po- 
tential in the radial direction is much weaker than that 
in the axial direction as uo r (t)/uo z <C 1, such that the 
motion of atoms in the z direction is essentially frozen 
to the ground state cp(z) of the axial harmonic trapping 
potential. 

Then we can separate the wave function as ^(r, t) = 
ijj(x,y,t)(p(z) to derive the 2D GP equation 

^ = -^^u J y4,+G V \^+ih{i ^, (i) 

with rj = J dz\(p(z)\ 4 / J dz\(p(z)\ 2 . Introducing the scales 
characterizing the trapping potential, the length, time, 
and wave function are scaled as 

x = a h x, y = a h y, t = i/u z , i/j $/a h y^Tra^, 

respectively, with ah = (h/mujz) 1 ^ 2 and ao is a constant 
length chosen to measure the s-wave scattering length. 



After the tilde is omitted, the 2D GP equation (pQ) is 
reduced to a dimensionless form as 

*W = ~\ V ^ + V( r >W + 9{r,t) \M 2 *P + in ^, (2) 

where the interaction strength g(r,t) = a(r,i)/ao, ft = 
VLq/uj z and the radial trapping potential can be written 
as V(r,i) = ou 2 (t)r 2 /2 with u{t) = uj r {t)/uj z . In what 
follows, we consider not only the harmonic potential like 
this but also an anharmonic potential. 

In order to find the exact vortex solutions to Eq. (j2]) 
with spatiotemporally modulated interaction, we first as- 
sume its exact solution as 

^(r, 0, t) = e iSe+i ^ r ^ p(r, t) U [R(r, t)] , (3) 

here S is the topological charge related to the angular mo- 
mentum of the condensate, p(r, t) denotes the amplitude 
of wave function and i?(r, t) is a intermediate variable 
reflecting the changes of main wave function U. Substi- 
tuting Eq. ((3]) into ((2]) and furthermore, letting U[R(r, t)] 
satisfy 

d 2 U/dR 2 + noli + piU 3 = 0, (4) 

where po and p\ are real constants, we can get 
a set of partial differential equations (PDEs) about 
p(r,t),R(r,t),</)(r,t),V(r,t) and#(r,£) as 

Rt + <j> r Rr = 0, 

2gp 2 ^p 1 R 2 = 0, 
pR r + 2 rp r R r + rp R rr = 0, 
rpcj) rr + 2 rcj) r p r + pcj) r + 2rp t = 0, (5) 

^-2(j) t -(j) r 2 +2ttS+-^-2V- p Rr 2 = 0, 

p r A pr 

here subscripts r and t mean the derivative of function 
with respect to r and t. If letting V(r,i) = co(t) 2 r 2 /2 — 
PqR 2 /2 (po = corresponds to harmonic potential) and 
0(r, t) = fir 2 + /2 (fi and fi are time-dependent func- 
tions, and fi is frequency chirp and fi is phase), then 
solving the set of PDEs (|5]) we get 

p(r,t) = e~ 2 t hdt e(re- 2 $ hdt ), (6) 

and 

re -2ffldt 

R(r,t)= / l/[B 2 (r)r]dr, (7) 

Jo 

and 

g(r,t) = -p 1 R 2 r /2p 2 J (8) 

where p\ is a parameter controlling the sign of in- 
teraction parameter g(r,t), 0(r) is defined by Whit- 
taker M and W functions [32|, i.e. 6(r) = 
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FIG. 1: (Color online) The spatiotemporally dependent in- 
teraction parameter g(r,t) in Eq. Q for parameter Q = 
0.7, Ai = 4,A 2 = 2,ci = c 2 = 1,5 = = 1000 and 

cj(t) = 0.028. (a) The radial structure of g(r,t) at t = 0. 
(b) The spatiotemporal structure of g(r,t). 



[c 1 M(A 1 /4A 2 ,5/2,A 2 r 2 ) + C2W(Ai/4A 2 , 5/2, A 2 r 2 )]/r, 
where Ai,A 2 ,ci,c 2 are nonzero constants and cic 2 > 0. 
In particular, the above f\ and / 2 satisfy the following 
two ordinary differential equations 



29. S- Aie" 4 S hdt - 2df 2 /dt = 0, 
uj 2 (t) + 4/ 2 + 2dfi/dt - \ 2 2 e- 8 $ hdt = 0. 



(9) 



When the parameter /xq = 0, the external potential 
is just harmonic form V(r,t) = uj 2 (t)r 2 /2, and we get 
explicit solution of Eq. (jl]) as !7(i?) = z/icn(z/i R + 
^o, V2/2)/ yfjli, where ^0 and ^1 are arbitrary constants 
and fi\ > 0. So the exact vortex solution to Eq. ((2]) is 



/Mi 



cn(z/ii? + zy , V2/2), (10) 



where functions p and R are given above. Here cn 
and sn (below) are Jacobi elliptic functions. When 
imposing the boundary conditions for vortex solution 
as lini| r |^o VK 7 "? <M) = li m |r|^oo VK r > ^M) = 0? we can 
get v = K(V 7 2/2),is 1 = 2nif(v / 2/2)/ J R(+oc,0), where 

K(s) = f* /2 [l - s 2 sin 2 r]-^ 2 dr is the first kind elliptic 
integral and n is a nonnegative integer. In this case, since 
111 should be positive and thus the interaction strength 
g(r,t) is negative corresponding to the condensates con- 
sisting of 85 Rb [lH or 7 Li atoms 34], [35j experimentally. 
The structures of the interaction parameter g(r,t) with 
respect to radial coordinate r and time t are demon- 
strated in Fig. [TJ It is observed that the interaction 
parameter is inverse Gaussian in r and periodic in t. 

When the parameter /iq 7^ 0, the external poten- 
tial becomes V(r,i) = u 2 (t)r 2 /2 — /ioR 2 /2 (an anhar- 
monic potential) as shown in Fig. [2j where there is 
a convex hull in the center of the harmonic poten- 
tial, and the anharmonic potential is periodic in time 
t. We get the exact solution of Eq. ((4]) as U(R) = 

y/2 (5 2 - /i )/ Misn (5R, y/fio/5 2 - l) , where /i /2 < 

5 2 < /Iq and pi < 0. So the exact vortex solution to 



FIG. 2: (Color online) The shapes of the anharmonic potential 
V(r,t) = uj 2 (t)r 2 /2 - fjL R%/2. (a) The radial structures at 
t — with parameters S — 1,5 = 7.4 (dashed line) and 
S — 5,(5 = 1172.4 (solid line), respectively, (b) The spacial 
structure at t = with parameter S = 1,5 = 7.4. (c) The 
spatiotemporal structure with parameter S = 1, 5 = 7.4. The 
other parameters are Q — 0.7, Ai — 4, A2 — 2, c\ — C2 — 
1,/zo = 62.2 and u(t) = 0.028. 



Eq. © is 
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l 2(5 2 -p ) 
jii 



pe< se+ ^ 2+ f*)sn (SR, VW<* 2 - l) , 

(11) 

where functions p and R are given above. When im- 
posing the boundary conditions for vortex solution as 
lini| r |^ VK^^M) = li m |r|^oo ^( r ? 0?^) = 0? we can get 
(5 = 2nK( y^/io/S 2 — l)/i?(oo,0). In this case, the pa- 
rameter /ii should be negative and thus the interaction 
strength g(r,t) is positive corresponding to the conden- 
sates consisting of 87 Rb [Hj or 23 Na atoms [I?} experi- 
mentally. 

Seen from the exact vortex solutions ([TO]) and ([TT]) . 
there exists two class of vortex states (distinguishing 
them with two quantum numbers which are radial node 
n and topological charge 5, also called angular momen- 
tum quantum number) corresponding to the harmonic 
potential (po = 0) and anharmonic potential (po 7^ 0) in 
attractive and repulsive BECs, respectively. In the fol- 
lowing, we will first examine the structures of these exact 
vortex solutions and then study the dynamic properties 
and stability of these vortex states under different situa- 
tions. 



STRUCTURES OF VORTEX STATES 

The structures of the exact vortex solutions ([TP]) and 
([TT]) can be controlled by modulating the frequency of 
the trapping potential and the spatiotemporal inhomo- 
geneous s-wave scattering length as seen from Eq. ([9]). 
Taking into account the feasibility of the experiment, we 
only consider the case of harmonic potential (po = 0) 
which corresponds to the attractive BEC as explained 
above. 

In real experiment, we assume an attractive 7 Li con- 
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FIG. 3: (Color online) The density distributions \^(x,y,0)\ 2 
(a)-(d), and the corresponding radial wave profiles (bottom) 
for the vortex solution ([TO]) of the attractive rotating BEC for 
topological charge 5 = 1 and various radial quantum num- 
bers. The parameters are Q = 0.7, jii = 1000, Ai = 4, A2 = 
2, ci = c 2 = 1, e = and u = 0.028. 



densate in the internal atomic state \F = l,rriF = 1) 
0, Ell trapped in an axis- symmetric disk-shaped po- 
tential, where the axial confinement energy Huj z is much 
larger than the radial confinement and interaction en- 
ergies, and the radial frequency of the trap is time- 
dependent which can be written as 



uj(t) = uj r (t)/uj z = ujq + ecos(cjit), 



(12) 



with < e < ujo. For e = 0, the radial frequency of the 
trap is time-independent and here, we choose the time- 
independent part of radial frequency uo r = (2tt) x 18 Hz 
and axial frequency uj z = (2tt) x 628 Hz as in [llj, so 
u(t) =ouo = 0.028. 

As shown in Fig. [3j we demonstrates the density distri- 
butions for different radial quantum number n with fixed 
topological charge S = 1 at t = 0, which is based on the 
exact vortex solution ([TP]) . The Fig. [3^a) corresponding 
to n = 1 is a lowest energy state and Figs. [3fb)-|3fd) 
corresponding to n = 2,3,4 are three excited states. In 
the Fig. [3fe), we show the radial wave profiles of the 
exact vortex solution ([TO]) at t = 0. It is clear to see that 
the number of ring structure of vortex solution increases 
by one with changing the radial quantum number n by 
one, which is similar to the quantum states of harmonic 
oscillator. 




Radius ( r ) 



FIG. 4: (Color online) The radial structures of the density 
distributions and phase diagrams for the vortex solution (|10j) 
of the attractive rotating BEC with radial quantum number 
n — 1 and different topological charges at t — 0. The insets 
are the corresponding phase diagrams and the other parame- 
ters are the same as Fig. 3. 
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FIG. 5: (Color online) Time evolution of the monopole mo- 
ment (r) in the attractive rotating BEC for different fre- 
quencies of the harmonic potential, (a) The radial quan- 
tum number is n = 1 and (b) n = 2. In both figures, the 
topological charge is S = 1 and the frequency of the trap 
is oj(t) = ojo + ecos(cjit) with parameters Q = 0.7, cjo = 
0.028, Ai = 4, A 2 = 2, m = 1000, a = c 2 = 1. 



One of the interesting properties for the exact vortex 
solution ([TO]) is shown in Fig. 2] by choosing different 
topological charge S with fixed radial quantum number 
n = 1. We can see that the density profiles of the vor- 
tex states become more and more localized with increas- 
ing the topological charge S due to the larger angular 
momentum for the higher topological charge S. More- 
over, vortex expands outwards with the increasing of the 
topological charges and so will obtain the larger angular 
momentum. 
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FIG. 6: (Color online) Time evolution of the density distribu- 
tions \ip(x,y,0)\ 2 and phase diagrams for the vortex solution 
(|1Q[) of the attractive rotating BEC for radial quantum num- 
ber n = 1. (a) Stable vortex for topological charge S — 1. 
(b)-(c) Unstable vortex for S — 2 and 3, respectively. For all 
cases, the domain is (x,y) £ [—15, 15] x [—15, 15]. The other 
parameters are Q = 0.7, Ai = 4, A2 = 2, /ii = 1000, c\ — 02 — 
l,e = and cjo = 0.028. 



Another interesting aspect of the condensate is to 
study the monopole moment [HI, [39| defined by (r) = 
J r|^| 2 dr which can be directly compared with experi- 
ments in BEC. In Fig. [5j we show the time evolution 
of the monopole moment for different oscillation frequen- 
cies (Ji, amplitude e in ([T2]h and different radial quantum 
numbers n = 1, 2 with fixed topological charge 5 = 1. 
It is seen that the monopole moment represents regular 
oscillation following the oscillating frequency of the trap 
when uji is large, but irregular oscillation when oji is small 
which can be observed experimentally. The amplitude 
and center position of oscillation of monopole moment 
for the case n = 2 are a little larger than that for the 
case n = 1, but their periods are the same. In particular, 
at e = i.e. time-independent trapping potential, the 
oscillation of the monopole moment is only determined 
by the spatiotemporally nonlinear interaction which also 
represents regular behavior in our studied case. 



STABILITY ANALYSIS 

It has been shown that attractive Bose condensates like 
85 Rb and 7 Li become mechanically unstable and collec- 
tively collapse 29|, [35| when the number of atoms in the 
condensate exceeds critical value N c . So it is important 
to produce the stable states in attractive Bose conden- 
sates. Saito and Ueda [3l| have demonstrated that a 
matter-wave bright soliton can be stabilized in 2D free 
space by causing the strength of interactions to oscillate 
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FIG. 7: (Color online) The dynamic instability of the vor- 
tex solution (|10[) of the attractive rotating BEC with time- 
dependent harmonic potential for topological charge S = 3 
and two different radial quantum numbers. The parameters 
are Q = 0.7, Ai = 10, A 2 = 5,/ii = 1000, a = c 2 = l,u>i = 
l,e = 0.02 and cj = 0.028. 



rapidly between repulsive and attractive by using, e.g., 
Feshbach resonance [29[. In previous work, we [23] have 
found an exact stable localized nonlinear matter wave 
in quasi-2D BEC with spatially modulated nonlinearity 
in harmonic potential. In this section, we investigate the 
dynamical stability of the exact vortex solutions ([TO]) and 
(pTj) by numerical simulation of Eq. (j2j). We show that 
only some types of the stable vortices (defined by radial 
and angular quantum numbers) can be supported by the 
spatiotemporally modulated interaction in this system. 

In order to elucidate the dynamical stability of the ex- 
act vortex solutions proposed in Section II, we conduct 
numerical experiments by solving Eq. ([2]) and take the 
exact vortex solutions ([TU]) and ([TT]) at t = as initial 
data. To begin with, we consider the attractive rotating 
BEC with harmonic potential at e = in ([T2]h which has 
exact vortex solution ([TO]) . In Fig. [6j we show the den- 
sity evolutions and phase diagrams of vortex solution ([TO]) 
as initial condition with radial quantum number n = 1 
and different topological charge S or angular momentum 
quantum numbers based on numerical simulation of Eq. 
([2]). It is shown that only when topological charge 5=1, 
vortex solution ([TO]) is stable against perturbation with 
an initial Gaussian noise of level 0.5%, but for topological 
charge S > 2, giant vortex solution (fT0|) will be unstable 
and split into single charge vortices and so destruct the 
ring structures. 

When the harmonic trap is time-dependent which cor- 
responds to e =^ and uj\ ^ in ([T2]h Fig. [7] shows 
the unstable dynamics and phase diagrams of the giant 
vortex solution (fTOj) with S = 3 and two different radial 
quantum numbers n = 1,2 for the attractive rotating 
BEC. It is observed that the time-dependent frequency 
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FIG. 8: (Color online) Time evolution of the density distri- 
butions \ip(x,y, 0)| 2 and phase diagrams for the vortex solu- 
tion ([11)) of the repulsive rotating BEC for topological charge 
S = 1 and radial quantum numbers n = 1 and 2, respectively, 
with an initial Gaussian noise of level 0.5%. Here the do- 
main is (x,y) £ [—25, 25] x [—25, 25], and the parameters are 
5 = 7.4 (top) and 6 = 14.9 (bottom). The other parameters 
are Q = 0.7, Ai = 4, A2 = 2, ci = C2 = 1, /ii = — 10, e = and 
cjo = 0.028. 
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FIG. 9: (Color online) Time evolution of the density distri- 
butions \ip(x, y, 0)| 2 and phase diagrams for the vortex solu- 
tion (|llj) of the repulsive rotating BEC for topological charge 
S — 5 and radial quantum numbers n — 1 and 2, respec- 
tively, with an initial Gaussian noise of level 0.5%. Here the 
domain is (x,y) G [—25,25] x [—25,25], and the parameters 
are 5 = 1172.4 (top) and 5 = 2344.9 (bottom). The other 
parameters are the same as Fig. 8. 



of trap affects the dynamics of the vortex significantly. 

Next we consider the repulsive rotating BEC in an- 
harmonic potential V{r,t) = ou 2 (t)r 2 /2 — hqR 2 /2 with 
/io 7^ shown in Fig. [2j which has exact vortex solution 
([TT]) . In Fig. [8]and[9j we demonstrate the density evolu- 
tions and phase diagrams of vortex solution ([TT]) as initial 
condition with different radial quantum numbers n = 1, 2 
and fixed angular momentum quantum numbers 5 = 1 
and 5 = 5, respectively. It is very interesting to note 
that when the radial quantum number n = 1, the exact 
vortex solution ([TT]) is always stable even for very large 
topological charge 5 = 5, which is very different from the 
attractive rotating BEC with harmonic potential where 
a stable region for vortex solution ([TO]) was found only 
for 5 = 1 as shown in Fig. El Our results have provided 
a very promising method for stabilizing the giant vor- 
tex having very large topological charge 5 > 2 which has 
been conjectured unstable |30[ by tuning the external po- 
tential and nonlinear interaction simultaneously in time. 
Numerical simulation shows that for the radial quantum 
number n > 1, the vortex solution ([TT]) is always unstable 
for any topological charge 5. 

Finally, we investigate the effect of the slightly asym- 
metrical potential to the stability of quantized vortices. 
To do so, we take the asymmetrical external trap as 
V(r, t) = u 2 (t) [(1 + e x )x 2 + (1 + e y )y 2 } /2 - fi R 2 /2, where 
parameters e x and e y describe small deviations of the 
trap from the axisymmetry. The ENS group [Io[ stirred 
a BEC of 87 Rb confined in this kind of magnetic trap us- 
ing a focused laser beam. Fig. 10 shows the evolutions of 



density profiles and phase diagrams of the vortex solution 
(11) at t = 1000 for the inhomogeneous repulsive rotating 
BEC in this slightly asymmetrical anharmonic potential 
with parameters e x = 0.02 and e y = 0.03. Here the 
quantum number n = 1 and topological charges 5 = 1 
and 5, respectively, and the level of the initial Gaussian 
noise is still 0.5%. It is seen that the quantized vortices 
with 5 = 1 and 5 in slightly asymmetrical anharmonic 
potential are still stable at t = 1000. 

We now give an experimental protocol to observe the 
above vortex states in future experiments. For the at- 
tractive interactions, we take 7 Li condensate in internal 



atomic state \F = l,mF = 1) [34|, |35|, containing about 
6.55 x 10 4 atoms, confined in a pancake-shaped trap with 
radial frequency uj r = (2ir) x 18 Hz and axial frequency 
uj z = (2tt) x 628 Hz [38|. Experimentally, this trap can 
be determined by combination of spectroscopic observa- 
tions, direct magnetic field measurement, and the ob- 
served spatial cylindrical symmetry of the trapped atom 
cloud [4l|. For the repulsive interactions, we take 87 Rb 
condensate in internal atomic state \F = l^mp = 1) 
(36| , containing about 8 x 10 5 atoms, confined in an an- 
harmonic potential which is a pancake-shaped harmonic 
trap with radial frequency uo r = (2tt)x 18 Hz and ax- 
ial frequency uj z = (2ir) x 628 Hz [38| plus a convex 
hull, see Fig. 2. The key step is how to realize the 
spatiotemporal variation of the scattering length. Near 
the Feshbach resonance [29|, the scattering length a s (B) 
varies dispersively as a function of magnetic field 5, i.e. 
a 3 {B) = a[l — A/(B — Bo)], with a being the asymptotic 
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FIG. 10: (Color online) Density profiles and phase diagrams 
for the vortex solution pip at t = 1000 for the repulsive 
inhomogeneous rotating BEC in slightly asymmetrical anhar- 
monic potential V(r,t) = 0J 2 (t)[(l + e x )x 2 + (1 + e y )y 2 }/2 - 



jjLoR r /2 with 



0.02 and 



0.03. Here the quantum 



numbers n = 1 and topological charges S = 1 and 5, respec- 
tively, and the level of the initial Gaussian noise is still 0.5%. 
The parameters are S = 7.4 (top) and S = 1172.4 (bottom). 
The other parameters are the same as Figs. 8 and 9. 



value of the scattering length far from the resonance, Bq 
being the resonant value of the magnetic field, and A 
being the width of the resonance at B = Bo. For the 
magnetic field in z direction with gradient a along x-y 
direction, we have B = [Bo + aBi(x, y, t)]e z . In this case, 
the scattering length is dependent on and time t. So 
in real experiments, we can use Feshbach resonance tech- 
nique to realize spatiotemporal variation of interaction 
parameters shown in Fig. 1. Finally, in order to observe 
the density distributions in Figs. 3 and 6-9 clearly in 
experiment, the atoms should be evaporatively cooled to 
low temperatures, say in the range of 10 to 100 nK. 



CONCLUSIONS 

In conclusion, we have investigated the quantized 
vortices in a rotating BEC with spatiotemporally 
modulated interaction in harmonic and anharmonic 
potentials, respectively. Two families of exact vortex 
solutions for the 2D GP equation are constructed 
explicitly by similarity transformation. It is interesting 
to see that a new series of stable giant vortex states with 
topological charge S > 2 can be supported by tuning the 
external potential and the spatiotemporally modulated 
interaction in this system. We hope that this paper will 
stimulate further research on quantized vortices and 
help to understand the behavior of nonlinear excitation 
in physical systems with spatiotemporally modulated 
interaction. 
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